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M3 1 (Andromeda) is our nearest neighbour spiral galaxy, and both its high mass and its low distance make it a source of 
interest for the indirect search for dark matter through y-ray detection. The ground based atmospheric Cherenkov telescope 
CELESTE observed M31 from 2001 to 2003, in the mostly unexplored energy range 50-500 GeV. These observations provide 
an upper limit on the flux above 50 GeV around 10~ 10 cm~ 2 s _I in the frame of supersymmetric dark matter, and more generally 
on any gamma emission from M3 1 . 
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Abstract. If dark matter is made of neutralinos, annihilation of such Majorana particles should produce high energy cosmic 
rays, especially in galaxy halo high density regions like galaxy centres. 

1. Introduction 

The presence of dark matter in the Universe has been known for decades. Since early measurements in galaxy clusters 
(Zwicky, I1937> . the mass distribution of the Universe has been studied at different scales with a focus on dynamical effects 
by means of galaxy star rotation curves, large-scale galaxy cluster dynamics (Ostriker et al.. U974> . Recent developments in ob- 
servational techniques in cosmology have resulted in independent estimates of the matter content of the Universe Q m /i 2 , through 
large-scale structure surveys (Hawkins et al.. l2003| : Loveday et al.. l2002t and measurements of the cosmic microwave back- 
ground radiation (the most recent being the WMAP mission (Spergel et al.. l2003» . Given standard cosmology, all suggest that 
most of the matter in the Universe is dark, cold and non-baryonic. Such hypotheses have led to the construction of the Cold Dark 
Matter (CDM) paradigm (Blumenthal et al., 1984 1: dark matter would be made of Weakly Interacting Massive Particles (WIMPs) 
which are neutral, stable and originating from the Big Bang era (Lee & Weinberg. 1 19771 . 

Supersymmetric theories (SUSY) (see for instance Nilles. [19841 offer an excellent WIMP candidate (Goldberg. 1 19831 . the 
neutralino, which is a mixture of the superpartners of the neutral Higgs bosons and of the electroweak gauge bosons. We will not 
discuss the case of extra-dimension phenomenological theories, which also provide interesting candidates (Servant & Tait. l2003l . 

The nature of dark matter is presently probed in both direct searches, by means of underground experiments that could detect 
elastic interactions of neutralinos with nuclei, and indirect searches using ground based or satellite telescopes to detect cosmic 
rays (gamma, leptons or hadrons) created by neutralino pair annihilations in galactic or extragalactic media (for review, see 
Bergstrom ROOOl . 

These different types of searches, together with collider experiments, are necessary to constrain in a wider view the quantum 
nature of dark matter, because they allow either additional or complementary surveys of the particle model parameter spaces. 

Searching for WIMP annihilation signatures with ground based y-ray telescopes leads to the question of the choice of targets. 
A good candidate will have a large amount of dark matter, and combine as big a density as possible, as small a distance from us 
as possible, and finally, will transit at high elevation in the experimental sky. The Galactic centre is a prime candidate except for 
being too close to the horizon for CELESTE, which is in the northern hemisphere. Instead, we have chosen M31 and the Draco 
dwarf galaxy for our searches. M31, which is the nearest giant spiral galaxy, is very massive (~ 1O 12 M ), and its star rotation 
curve indicates a large amount of dark matter. Draco, a neighbour dwarf spheroidal galaxy dominated by a dark component 
(Kleyna et al.. l2001> . is also a very good candidate but our attempts to study it were foiled by bad weather. 

In this paper, we present the result of searches for y-ray emission from M31 with the CELESTE telescope. In section |2] we 
review the predictions made for the observations (more details in Falvard 2004, hereafter F04), for which we considered CDM in 
the frame of minimal supergravity (mSUGRA) phenomenology. The impact of the astrophysical modelling is briefly revisited, as 
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well as possible consequences of non-standard cosmologies. Focusing on the experimental techniques, we present in section[3]the 
method we use to search for a y-ray signal, with explicit comparisons to Crab data. In the absence of a detection, 2<x confidence 
level limits are computed for all studied SUSY models. 



2. Gamma-ray flux predictions for supersymmetric annihilating dark matter in M31 

Under the assumption that neutralinos have an isotropic and homogeneous velocity distribution, which is likely to be the case in 
halo centres where the WIMP velocity is expected to be low, the averaged y-ray flux due to their annihilation, integrated above 
an energy threshold and within the solid angle AO. can be written as: 



1 N y (E th ) <crv> f f 1 N 7 (E & ) < crv > 

<D(£th) = ~. P ( s ) ds d ^ = ~ A 7T-2 C 1 ) 



We thus decouple astrophysics modelling from SUSY contributions. The first part of the right hand term is related to particle 
physics, via the thermally averaged product of the cross section <x with the velocity v producing N y (E t h) photons of energy 
E > E t h, and the neutralino mass m Xo . The second term refers to the (squared) halo density profile p integrated within an 
experimental field of view of angular radius 9 along the line of sight ds. At the same time, we define 1,(0) as the astrophysical 
factor of the flux. 

As y-rays result mainly from hadronization of annihilation final states (mainly quarks and gauge bosons), their spectral 
shape mainly takes its origin in the decay of n°. It has been shown by several authors (see for instance Bergstrom et al.. U998| 
or Tasitsiomi et al., 2002) that such a spectrum can be fitted or modelled with respect to the neutralino mass. Therefore, the 
differential spectrum above a threshold energy E± can be written as follows: 



d<S> 
dE 



—(E > E th ) = (D(£ th ) x f(E, m xa ) (2) 



where f(E,m Xo ) is the spectral shape derived from the SUSY model and normalized such that f(E,m Xa )dE — 1, so that 
<5(£' t h) is the integrated spectrum above an energy threshold E^. This expression will be useful when we interpret the M31 data 
collected by CELESTE. 



2.1 . Halo modelling 

M3 1 is a late-type Sb spiral galaxy, which lies at a distance of about 675 kpc, and is observable from the Northern hemisphere (RA 
= 10.68°, DEC = 41.27°). A study by Braun dl991l l. based upon the analysis of HI data and a model-independent reconstruction 
of the velocity field, showed that the star rotation curve arises naturally by considering two optically traced mass components: 
a bulge, with a total mass of (7.8 ± 0.5) x 1O IO M , and a disk of (1.22 + 0.05) x 1O U M within 28 kpc. Nevertheless, it seems 
that the star mass-to-light ratios used in this paper, T 'bulge =6.5 and T^isk = 6.4 (solar units in blue band), are over-estimated (cf. 
F04). 

By lowering the bulge and disk contributions, that is Tbuige — 3.7 and Tdisk - 2.5 as indicated by F04, we assumed that a dark 
halo significantly accounts for the gravitational potential. Let us consider the following CDM density profile: 

PCDMM=Po(^) y (^) , (3) 

where r ( ) usually stands for a core radius and a a scale length. An Navarro-Frenk- White profile (Navarro et al., |T996| l, i.e. with 
y = 1, a = 1 and e — 2, fits the rotation curve, and has its parameters entirely determined by the previous mass-to-light ratios. 
The fitted values are po ~ 0.07GeV.cnT 3 , ro ~ 20kpc and a ~ 5kpc. The resulting contribution to the y-flux within CELESTE'S 
field of view (9 = 5mrad ~ 0.3°, corresponding to a 3.5 kpc radius) has been calculated: 

I(6» = 5mrad) = 3xl0 19 GeV 2 cirT 5 . (4) 

This result depends strongly on the central tail of the dark halo, but is rather conservative since calculated with a profile. 



2.2. Probing the SUSY parameter space 

We choose the minimal supergravity framework (mSUGRA) to scan over the SUSY parameter space. In this frame, a SUSY 
model can be defined at the unification scale with 5 parameters: the unified scalar mass mo, the unified gaugino mass nt\/2, the 
Higgs vacuum expected value ratio tan(/?), the unified trilinear coupling Aq and the sign of the mixing parameter of the Higgs 
superfields yu. We use an interface between the public codes Suspect (Djouadi et al.. l2002> and DarkSusy (Gondolo et al.. l2004> 
to compute SUSY masses, annihilation rates and relic densities for various random models. The constraints on these models 
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Fig. 1. Integrated expected y-flux from M31 as a function of the energy threshold, for two selections of low (100-1 10 GeV) and 
high (500-600 GeV) neutralino masses (Q Xo h 2 e [0.086,0.14]), showing the mass-dependence of the flux, and the complemen- 
tarity between EGRET and CELESTE. 



come from standard accelerator limits, and we select a rather large range for the relic density (Cl Xo h 2 e [0.05, 0.14]) 1 , according 
to the WMAP measurement Cl xo h 2 = 0.113 ± 0.009 (Bennett et al.. 120031 . 

Combining the resulting y-spectrum with the astrophysical factor E, we calculate the integrated flux as a function of the 
threshold energy, as shown in Fig.^ m which we have selected only two groups of WMAP compatible neutralinos (respectively 
small and high masses) in order to exhibit the mass-dependence of the expected flux. The SUSY models plotted there are 
characterized by a high value of tan/3 (typically > 30), for which the production of y-rays is more efficient (due to high branching 
ratio in b quarks). On the same plot, we represent the upper limit provided by the EGRET collaboration (Blom et al.. U999t and 
CELESTE'S sensitivity, ~ 6 x 10 _11 ph.crrr 2 s _1 , estimated for 50 hours of observation. This figure illustrates the mass-and/or- 
energy experimental complementarities. 

2.3. Other contributions 

Although the predicted fluxes are low (~ 10~ 13 crrr 2 s~' at a threshold of ~ 50 GeV), ~ 3 orders of magnitude smaller than for 
the Crab nebula (therefore far from CELESTE'S sensitivity), several effects could enhance them. 

First of all, dark matter substructures, the so-called clumps, arise naturally in the hierarchical scheme of galaxy formation, 
and simulations of the non-linear regime of collapse allow a semi-theoretical study of their statistics and structure (Moore et 
al.. U999> . Such local overdensities should induce extra neutralino annihilations, and translate to an additional factor to the flux. 
Although it is rather difficult to estimate how clumpy the galaxies remain today, this enhancement factor was until recently 
supposed to be smaller than ~ 10 (Stoehr et al.. l2003t . Nevertheless, a recent study by Diemand et al. ( 120051 suggests that about 
50% of a Milky- Way-like galaxy mass is possibly bound to dark matter substructures, whose mass range spreads from 10~ 6 up 
to 10 7 M o . The authors claim that about ~ 10 15 of such substructures may have survived against gravitational disruption, leading 
to a boost factor of over two orders of magnitude compared to the smooth contribution. M3 1 being very similar to our galaxy, 
this statement should also stand for that source. 

As another possible astrophysical effect, the supermassive black hole at the centre of M31 could raise the central halo 
profile up due to adiabatic accretion (Gondolo & Silk, 1999 1. There are other interesting mechanisms involving baryons to 
enhance the dark matter density. Gnedin et al. J2004> sketch such responses of halos to condensation of baryons, while Bertone 



1 The upper limit is given by the WMAP result plus three sigma. Higher values for the relic density are not that interesting because they 
correspond to lower values of the annihilation cross-section. 
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et al. J2005> propose a possible enhancement of the WIMP annihilation rate due to the presence of intermediate-mass black holes. 

Beside those astrophysical effects, some recent developments in the frame of theoretical cosmology have focused on the 
quintessence scheme (Caldwell et al.,[l998 1 to solve the so-called coincidence problem (the fact that Qa ~ ^matter today). Such a 
quintessential field could undergo a kination regime in the early universe (Salati. l2"003> . so that its kinetic energy dominates over 
its potential. In this regime, the expansion rate of the universe is enhanced and the thermal history of neutralinos is consequently 
modified: the decoupling of neutralinos can take place more rapidly at earlier times. Therefore, the WMAP constraint leads 
to a higher neutralino annihilation cross-section. This means that this phenomenon rehabilitates SUSY models for which relic 
densities are too low, when calculated in standard cosmology. Salati J2003> shows that the relic density enhancement can be 
parameterized by: 



?7o = po,o/Py,o, where the 0-index refers to a temperature of 1 MeV, and p<p (respectively p y ) is the quintessence (photon) energy 
density. The upper limit on 770 comes from Big Bang Nucleosynthesis stages that should not be perturbed by the kination regime 
(Yahiroet al.. !2UU2l . 

According to this cosmology, higher annihilation rate models are required, which therefore means that the y-ray production 
is enhanced. This effect appears in the final results showed in Fig. [7] for a small sample of selected SUSY models. 

Finally, we emphasize another interesting effect coming from Affleck-Dine baryogenesis in SUSY (Fujii et al.. l2004> . which 
yields natural matter-antimatter asymmetry in the early universe. In such a scenario, meta-stable particles result from oscillations 
in flat directions of the scalar potential, carrying baryon and/or lepton number, namely Q-balls. These Q-balls can have a 
lifetime long enough to decay after the freeze-out of neutralinos. This induces a non-thermal production of neutralinos, and thus 
enhances their relic density. This also requires, as previously, higher neutralino annihilation rates to not overdose the universe. 

Therefore, although standard conservative predictions are not that optimistic, all these putative contributions may increase 
the y-flux from M3 1 significantly. This further motivates observations of such a source with CELESTE, keeping in mind that 
CDM could be something besides SUSY. 

3. Observations of M31 with CELESTE 

3.1. The CELESTE experiment 

CELESTE (Pare et al. I2002I de Naurois et al. I2002t is an atmospheric Cherenkov telescope detecting y-rays above ~ 50 GeV 
(the experiment shut down in June 2004). Reaching such a low energy threshold required a large light collection area, achieved 
by exploiting the mirrors of a solar plant. These mirrors are used to sample the arrival time and photon flux of the Cherenkov 
wavefront generated by atmospheric showers initiated by cosmic rays at many points in the light pool. In contrast to the imag- 
ing technique (Weekes. ll988> the sampling technique uses information on the shape of the wavefront for hadron rejection, as 
described below. 

The CELESTE experiment uses 53 heliostats (40 until 2001) of the Themis former solar plant (French Pyrenees, 42.50°N, 
1.97°E, altitude 1650 m). Each heliostat (54 m 2 ) reflects the light onto the secondary optics, located at the top of a 100 me- 
ter tower, focussing the light onto a single photomultiplier (PMT) for each heliostat. The PMT signals are sent to the trigger 
electronics and to the data acquisition system where they are digitized by ~ 1 GHz flash analog-to-digital converters (FADCs). 

The mean altitude of the maximum Cherenkov emission for y-ray induced showers is around 1 1 km above the site. The 
heliostats are aimed at this altitude in the direction of the source under study to enhance light collection. The observations are 
made in the On-Off tracking mode: the observation of the source (On) is followed or preceded by an observation at the same 
declination offset in right ascension by 20 minutes. The latter is used as a reference for the cosmic-ray background and the signal 
is given by the difference between On and Off, after analysis cuts. 

M31 is a special source for CELESTE in that its blue magnitude is about 4.3, and ~ 5.4 if integrated in the +5 mrad field 
of view of CELESTE (de Vaucouleurs, [1 95 8 ) . Hence, pointing On-source increases the PMT illumination compared to the Off- 
source data. The same problem arises to a lesser degree for the study of the blazar Mrk421, due to the presence of the mg = 6.16 
star 5 1 U Ma in the same field-of-view. On-off illumination differences introduce biases at the trigger level and in the recorded 
data which fake a signal if not handled properly. We remove these biases by "padding" the Off-source data with extra background 
light as part of the data analysis, and by making a pulseheight cut 10% above the hardware trigger threshold, as described in (de 
Naurois et al.. l2002> and updated in (Manseri. l2004> . 

In the following, all Monte Carlo simulations are performed at the transit position of sources in the Themis sky, unless 
specified. We have used the atmospheric shower simulator Corsika (Heck et al.. U998> for our Monte Carlo studies. Moreover, 
stellar photometry studies using the PMT anode currents provided an improved description of our optics in the detector simulation 
(Smith & Brion. l2004l 




(5) 
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Fig. 2. Left: normalized distributions of the discriminating analysis variable £ for a simulation of y-rays (solid line) from M3 1 
following a CDM-like spectrum (m Xo = 500 GeV), with respect to M31 Off-data (2001-03, black markers). Right: number 
distribution of £ for a E^-spectrum simulated at the Crab transit (solid line), with respect to the On-Off difference for a data 
sample of the Crab nebula (taken between 2002 and 2004, black markers with error bars). 

3.2. Hadron rejection 

In our energy range, the Cherenkov wavefront from y-ray induced showers is, on average, more spherical than for showers 
initiated by hadrons (that is, by charged cosmic rays, mainly protons). The CELESTE field of view is small compared to the 
angular extent of the showers, which lessens the difference, but efficient hadron rejection is still possible. 

Just above the trigger threshold, the Cherenkov signal for many heliostats is comparable to the fluctuations of the night sky 
background light, so we sum all the signals instead of using each individually. Summing the signals implies compensating for 
the propagation delays, which requires knowledge of the shower core position when assuming a spherical Cherenkov wavefront. 
As described in Manseri 2004 this is done by scanning the plane at 1 1 km and evaluating the ratio H/W for each position of the 
scan, where H and W are the height and the width of the summed signal. The largest value of the ratio H/W yields our measure 
of the shower core position. The sphericity of the wavefront is estimated by how much the ratio H/W decreases when estimated 
200 m away from the shower core position. This relative decrease, called is shown in Fig.|2j-left for an Off observation and 
for a simulation of a y-ray spectrum. As expected, because of their Cherenkov wavefront sphericity, y-ray showers have lower £ 
values than hadronic showers. 

The On-Off difference of this relative decrease is shown in Fig. [JJright for a sample of Crab nebula data, which was taken 
between 2002 and 2004 with the same experimental setup as for M3 1 . It exhibits a clear excess due to the y-ray emission from 
the Crab nebula. Requiring £ < 0.35 leads to a 13.5<x excess and a sensitivity of 6.5a/ yfh. 

3.3. Observations, data selection and signal searches 

M31 was observed with CELESTE from 2001 to 2003, and 68 On-Off pairs were collected (~22 hours of On-source data). 
Nevertheless, variations in atmospheric conditions are known to cause systematic shifts in the On-Off difference, so we applied 
a selection based on criteria requiring stable detector operation (characterized by PMT anodic current and trigger rate stability). 
This selection reduced the data set to 6.5 hours because of bad weather conditions at Themis since 2001. 

The On-Off difference of £ for the M31 data is shown in Fig.|5J-left. No evidence of an excess can be found, and the On-Off 
difference is -0.75<x when requiring f < 0.35. 

SUSY annihilating dark matter could produce soft spectra or peaked signals, so we search for a signal in various energy 
bands. Our energy reconstruction uses y-ray simulations with fixed energies, and is based on the correlation between the total 
charge recorded with the FADCs and the true (i.e. Monte Carlo) energy, for a given impact parameter. (The shower core position 
at 1 1 km, obtained by maximizing H/W as described above, gives the impact parameter on the ground assuming that the gamma 
ray comes from the source under study.) Figure |4] (left) shows the mean charge per heliostat versus the reconstructed impact 
parameter, for different y-ray energies. We use the charge distributions at each energy and impact parameter to build the inverse 
function, that is, the function predicting the energy from the observed charge and the reconstructed impact parameter in the range 
of [0 - 120] m. To check the function, we inverse it again and compare it to the original points, also shown in Fig.0](left). The 
right panel of Fig. 0]shows the energy bias and resolution curves obtained from this estimation method. The energy resolution is 
below 30% with a small bias of ~ 5%, adequate to search for an excess in the data. 

The On-Off difference of the measured energy distribution, after analysis cuts, is shown in Fig. [3jright, in which the 50 
GeV binning covers at least one standard deviation of the energy resolution function. No excess has been found. All results are 
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£ Estimated energy E 7 (GeV) 

Fig.3. Left: On-Off distribution of £ for M31 data (2001-03, black markers with error bars), with respect to the same simulated 
spectrum as Fig.|3-left (solid line). Right: On-Off distribution of the estimated energy for all events within analysis cuts. 




Impact parameter (m) True energy (GeV) 

Fig. 4. Left: The crosses are the mean charge per heliostat (arbitrary units) for y-rays simulated with energies of 50, 100, 150, 
200, 300 and 400 GeV (bottom-top), versus the reconstructed core position (i.e. the impact parameter). The solid curves show 
the charge obtained from the functions used to predict the energy. Right: resolution (solid upper curve) and bias (lower dashed 
curve) for the energy reconstruction method with respect to the true energy (same Monte-Carlo set-up as in plot at left). 



Analysis level 


On-source evts 


Off-source evts 


On-Off 


Significance (No-) 


Raw data 


463520 


462327 


1193 


-0.59 


Analysis cuts 


10615 


10740 


-125 


-0.75 


E meas < 100 GeV 


6101 


6167 


-66 


-0.53 


100 < £ raclls < 200 GeV 


3143 


3197 


-54 


-0.61 


200 < £ mcas < 300 GeV 


870 


824 


46 


1.02 


300 < £ mcas < 400 GeV 


260 


269 


-9 


-0.35 



Table 1. Final statistics resulting from M31 data analysis. The analysis cuts are the following: we impose a software trigger 10% 
higher than the hardware one, and we apply a cut £ < 0.35 (see text for details) ; for signal searches in energy bins, we add a 
selection on the reconstructed impact parameter - < 120 m - according to the validity range of our energy reconstruction method. 
No significant excess appears over the whole sample, nor within different ranges of energy. 



summarized in Tabled in which a search for an excess within bins of 100 GeV, much larger at low energy than the resolution, is 
presented. 

3.4. Stability and upper limit 

We have also studied the stability of this result, as a check of possible systematic effects. We show in Fig.[5]how the significance 
remains stable when varying the cut value on the main discriminating variable ^. As a comparison, the same exercise has been 
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Fig. 5. Significance level according to varying the cut on £ (after raising the threshold by 10%). This figure shows the study of 
stability of the analysis as well for Crab as for M3 1 data. A signal evidence is clearly seen for the Crab, whereas the significance 
level remains fiat for M3 1 . 



done for a Crab dataset (same figure), and the y-ray signal clearly exhibits a bump in significance up to No- ~ 13.5 around 
£ < 0.35. 

On-Off difference measures the integral over energy of the theoretical spectrum of equation [2] convoluted with the effective 
area, times the observation time. As the statistics collected is compatible with the absence of a signal, we have only measured the 
background and its fluctuations in direction of M3 1 . This can be translated to an upper limit on a flux coming from that source, 
given a normalized theoretical spectral shape /(£). For a No- upper limit, and given an experimental energy threshold E^, any 
integrated flux above this energy should be bounded like: 

(D(£ > E th ) < N<r J Nbkgd (in cnrV 1 ), (6) 

T obs J E h Jl(E)f(E)dE 

where SNbkgd — -J2 X Noff is the measured background RMS, T„i, s is the total exposure time and 31{E) stands for the energy- 
dependent effective detection area. The latter is determined by means of simulations, and is plotted in Fig. |6jleft. According to 
equation|21 note that the spectrum f(E) is mass-dependent in case of neutralino annihilation, so that the limit should depend on 



the mass. For the sake of simplicity, we will use the parametrization given in Tasitsiomi et al. (2002 1 for the spectral shape, that 
is: 

6V^h /10 5 3x 2 +6x-r 



fiE) 5m w (V^-D 4 U 12 *3/2 j ( ? ) 

where x = E/m Xo and x t h = E±/m Xo . This spectrum depends on the mass, but assumes that y-rays come only from n° decay. 
Therefore, it does not take into account all specificities of SUSY models. Nevertheless, it is sufficient for our purpose. This 
parametrization is illustrated in Fig.[6j-right for two masses, 250 and 500 GeV, and compared with a l/E 2 power law spectrum. 

The threshold is set to 50 GeV, taken from the effective area shown in Fig. [6] (left). Given this threshold, a flux limit can 
be computed for each neutralino mass, using equation[6] The result is presented in Fig.0 where the predicted integrated fluxes 
of y-rays above 50 GeV are plotted with respect to neutralino masses. The averaged limit in a mass range of [50 - 700] GeV 
lies around 10 _10 ph.cm. _2 s _1 . We emphasize that this is the first experimental result in the energy range 50-500 GeV, and is 
complementary to those provided by EGRET (Blom et al., 1999) and HEGRA (Aharonian et al.. l2003> observations of M31. 



4. Discussion and conclusion 

Observations of M3 1 with CELESTE provide a 2cr upper limit on the y-flux above 50 GeV, depending on the expected spectrum. 
This limit, around 10 _10 ph.cm _2 s _1 , is quite far from the SUSY parameter space, but significantly constrains combinations of 
different enhancement factors discussed in § 12.31 (which are also likely to be excluded by EGRET limits, depending on the 
neutralino mass), and also any other model of annihilating dark matter besides SUSY. Whereas these observations have been 
motivated by indirect searches for SUSY CDM, this result yields a general astrophysical result: the first observation of a spiral 
galaxy in this energy range, somehow constraining y-ray emission from this class of objects. 

Any gamma ray detection from a galaxy like M3 1 would be difficult to interpret in terms of dark matter annihilation. Spiral 
galaxies are known sites of non-thermal processes and cosmic ray acceleration, and the relevant physical mechanisms are not yet 
well understood. In this sense, a Dwarf Spheroidal galaxy like Draco is a very promising source for indirect detection, given it is 
clearly dominated by the dark matter component. Unfortunately, we have too few data on Draco to perform a relevant analysis. 
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Fig. 6. Left: effective detection area - dashed line for trigger, full after event selection - as a function of the true energy, based on 
simulations of y-rays coming from the transit position of M3 1 . Right: normalized spectra between 50 and 500 GeV for a power 
law of index -2, and for formula of equation^with m x o = 500 (dashed curve) and 250 GeV (dotted curve). 

However, the upcoming generation of y-ray instruments will undoubtedly further constrain dark matter models and halo 
models for various astrophysical sources. These searches are not only complementary to future particle physics experiments, but 
also very important to understand how the question of dark matter is connected to the particle content of the Universe. 
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